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We investigate hyperfine induced electron spin and entanglement dynamics in a system of two 
quantum dot spin qubits. We focus on the situation of zero external magnetic field and concentrate 
on approximation-free theoretical methods. We give an exact solution of the model for homogeneous 
hyperfine coupling constants (with all coupling coefficients being equal) and varying exchange cou- 
pling, and we derive the dynamics therefrom. After describing and explaining the basic dynamical 
properties, the decoherence time is calculated from the results of a detailed investigation of the short 
time electron spin dynamics. The result turns out to be in good agreement with experimental data. 
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I. INTRODUCTION 

Quantum dot spin qubits are among the most promis- 
ing and most intensively investigated building blocks 
of possible future solid state quantum computation 
systemsiiS. One of the major limitations of the deco- 
herence time of the confined electron spin is its interac- 
tion with surrounding nuclear spins by means of hyper- 
fine interaction-"^-. For reviews the reader is referred to 
Refs.^^"^^. Apart from this adverse aspect, hyperfine in- 
teraction can act as a resource of quantum information 
processingi^"— . For the above reasons it is of key interest 
to understand the hyperfine induced spin dynamics. 

Most of the work into this direction, for single as well as 
double quantum dots, has been carried out under the as- 
sumption of a strong magnetic field coupled to the central 
spin system. This allows for a perturbative treatment or 
a complete neglect of the electron-nuclear "flip-flop" part 
of the Hamiltonian, yielding great simplification'^'ii^"— . 
In the present paper we consider the case of zero mag- 
netic field where such approximations fail, and we there- 
fore concentrate on exact methods. 

In the case of a single quantum dot spin qubit the 
usual Hamiltonian describing hyperfine interaction with 
surrounding nuclei is integrable by means of Bethe ansatz 
as devised by Gaudin several decades ago^"—. In the 
following we shall refer to that sytem also as the Gaudin 
model. Nevertheless exact results are rare also here be- 
cause the Bethe ansatz equations are very hard to han- 
dle. Hence there are mainly three different routes in or- 
der to gain some exact results: (i) Restriction of the ini- 
tial state to the one magnon sector— li., (ii) restriction to 
small system sizes enabling progress via exact numerical 
diagonalizationsiii^S, and (iii) restrictions to the hyper- 
fine coupling constants- 1''^-. In the present paper we will 
follow the third route and study in detail the electron spin 
as well as the entanglement dynamics in a double quan- 
tum dot model with partially homogeneous couplings: 
The hyperfine coupling constants are chosen to be equal 
to each other, whereas the exchange coupling is arbitrary. 
Although the assumption of homogeneous hyperfine con- 



stants (being the same for each spin in the nuclear bath) 
is certainly a great simplification of the true physical sit- 
uation, models of this type offer the opportunity to ob- 
tain exact, approximation-free results which are scarce 
otherwise. Moreover, such models have been the basis 
of several recent theoretical studies leading to concrete 
predictionsiir— 

The paper is organized as follows: In Sec. |TT] we in- 
troduce the Hamiltonian of the hyperfine interaction and 
derive the spin and entanglement dynamics for homo- 
geneous hyperfine coupling constants. In Sec. IIIII we 
study the spin and entanglement dynamics for different 
exchange couplings and bath polarizations. For the com- 
pletely homogeneous case of the exchange coupling being 
the same as the hyperfine couplings we find an empirical 
rule describing the transition from low polarization dy- 
namics to high polarization dynamics. The latter shows 
a jump in the amplitude when varying the exchange cou- 
pling away from complete homogenity. This effect as well 
as features like the periodicity of the dynamics are ex- 
plained by analyzing the level spacings and their con- 
tributions to the dynamics. In Sec. IIVI we extract the 
decoherence time from the dynamics by investigating the 
scaling behaviour of the short time electron spin dynam- 
ics. The result turns out to be in good agreement with 
experimental findings. 



II. MODEL AND FORMALISM 

The hyperfine interaction in a system of two quantum 
dot spin qubits is described by the Hamiltonian 
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where Jex denotes the exchange coupling between the 
two electron spins Si, S2, and Aj, are the couphng 
parameters for their hyperfine interaction with the sur- 
rounding nuclear spins Ii. In a realistic quantum dot 
these quantities are proportional to the square modulus 
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FIG. 1: (Color online) Spin dynamics for |qi) = \-fii\. 
),|r+),|ro) and an even number of spins. The number of 
down spins in the bath is Nd = 20 in all plots, yielding polar- 
izations pb ~ 5% — 30%. Note that the time unit is rescaled ac- 
cording to the number of bath spins. We see periodicity with 
TT. For = ITb) and = 58 we count the number of local 
extrema on one period and find N—2ND + i ~ 58 — 40+1 = 19 
as expected. 



of the electronic wave function at the sites of the nuclei 
and therefore clearly spatially dependent 

Al^Av\2p^n)\\ (2) 

where v is the volume of the unit cell containing one nu- 
clear spin and (fi ) is the electronic wave function of 
electron j = 1,2 at the site of i-th nucleus. The quan- 
tity Ai denotes the hyperfine coupling strength which 
depends on the respective nuclear species through the 
nuclear gyromagnetic ratio It should be stressed that 
these can have different lengths. In a GaAs quantum dot 
for example all Ga and As isotopes carry the same nu- 
clear spin li = 3/2, whereas in an In As quantum dot the 
In isotopes carry a nuclear spin of li = 9/3ii. In any 
case the Hamiltonian obviously conserves the total spin 
J = S + T, where S = Si + S2 a.nd I = J2f=i 

The model to be studied in this paper now results by 
neglecting the spatial variation of the hyperfine coupling 
constants and choosing them to be equal to each other 
A} = = A/N . Variation of the exchange coupling be- 
tween the two central spins J^x then gives rise to an in- 
homogeneity in the system. Hence the two electron spins 
are interacting with a common nuclear spin bath. More- 
over, if small variations of the coupling constants would 
be included, degenerate energy levels would slightly split 
and give rise to a modified long-time behavior of the sys- 
tem. In our quantitative studies to be reported on below, 
however, we focus on the short-time properties where de- 
coherence phenomena take place. Indeed, in section IIVI 
we obtain realistic T2 decoherence time scales in an al- 
most analytical fashion. In consistency with the homoge- 
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FIG. 2: (Color online) Spin dynamics for |ai) = \-\[\!y 
) , \T+), \To) and an odd number of spins. The number of down 
spins in the bath is Nn = 20 in all plots, giving polarizations 
Pb ~ 2% — 30%. In contrast to the case of an even number of 
spins we see periodicity with 27r. For \ai} — IffJI) and A*' = 45 
we count the number of local extrema on half the period and 
find N - 2Nd + 1 = 45 - 40 + 1 = 6 as expected. 



nous couplings we choose the length of the bath spins to 
be equal to each other. For simplicity we restrict the 
nuclear spins to = 1/2. We expect our results to 
be of quite general nature not strongly depeding on this 
choice^*. Note that both, the square S"^ of the total cen- 
tral spin as well as the square P of the total bath spin 
are separately conserved quantities. 

Considering the two electrons to interact with a com- 
mon nuclear spin bath as in our model corresponds to a 
physical situation where the electrons are comparatively 
near to each other. This leads to the question whether 
our model is also adapted to the case of two electrons 
in one quantum dot, rather than in two nearby quantum 
dots. Assuming perfect confinement, in the former case 
one of the two electrons would be forced into the first 
excited state, which typically has a zero around the dot 
center. Thus, the coupling constants near the very center 
of the dot would clearly be different for the two electrons. 
Therefore our model is more suitable for the description 
of two electrons in two nearby quantum dots than for the 
case of two electrons in one dot. 

Let us now turn to the exact solution of our homoge- 
neous coupling model and calculate the spin and entan- 
glement dynamics from the eigensystem. In what follows 
we shall work in subspaces of a fixed eigenvalue of J^. 
Thus, the expectation values of the x- and j/-components 
of the central and nuclear spins vanish, and we only have 
to consider their z-components. 

If all hyperfine couplings are equal to each other A] = 
Af = A/N, the Hamiltonian ([T]) can be rewritten in the 
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FIG. 3: (Color online) Spin dynamics for |qi) — and 
Nn = 20, resulting in pt ~ 6% — 30%. If Jex is an odd 
multiple of A/2N we see periodicity with 2n. 



FIG. 4: (Color online) Spin dynamics for joi) = 
(1/^) (21-^4) + 314-^)) and No = 20, resulting in pt ^ 
6% - 30%. 



following way 

H — Hhom + ^Jex ~ -Si • 52 (3) 

with 

i/ho.„ -^^ij"- Sf -Sl^ . (4) 

Omitting the quantum numbers corresponding to a cer- 
tain Clebsch-Gordan decomposition of the bath, the 
eigenstates are labelled by J, m, S associated with the 
operators J^,J^,S'^. The two central spins couple to 
S = 0,1. Hence the eigenstates of H are given by triplet 
states \J,m, 1), corresponding to the coupling of a spin 
of length one to an arbitrary spin, and a singlet state 
I J, m, 0). The explicit expressions are given by (|AH IA3| 
IA4[) in appendix A. 

The corresponding eigenvalues read as follows: 

H\I+l,m,l) = (^^i+i^yi+i^rn,l) (5a) 

H\I,m,l} = (^-4) l^''"'!) (5b) 

H\I~l,m,l) = 1,^,1) (5c) 

H\I,m,0) = -^J,^\I,m,0) (5d) 

Now we are ready to evaluate the time evolution of the 
central spins and their entanglement from the eigensys- 
tem of the Hamiltonian. We consider initial states \a) of 
the form |a) = |ai)|Q;2), where \ai) is an arbitrary cen- 
tral spin state and \a2) is a product of N states |t), 
The physical significance of this choice becomes clear 
by rewriting the electron-nuclear coupling parts of the 



Hamiltonian in terms of creation and annihilation oper- 
ators: 

S^Ij = I {S+Ii + Sri+) + Stl^ (6) 

Obviously the second term does not contribute to the dy- 
namics for initial states which are simple product states. 
Hence by considering initial states of the above form, we 
mainly study the influence of the flip-flop part on the dy- 
namics of the system. This is exactly the part which is 
eliminated by considering a strong magnetic field like in 
Refs.^^S-^. 

As the 2^ dimensional bath Hilbert space is spannend 
by the P' eigenstates, every product state can be written 
in terms of these eigenstates. If No < N/2 is the number 
of down spins in the bath, it follows 

^ ^ ,„ , TV TV 
\l...lt...^)^Y.Y.ci''^\--k,--Nn,{S^}), (7) 

where the quantum numbers {Si} are due to a certain 
Clebsch-Gordan decomposition of the bath. In ([7]) we 
assumed the first A'^^i spins to be fiipped, which is no loss 
of generality due to the homogeneity of the couplings. For 
the following discussions it is convenient to introduce the 
bath polarization pi, = (TV — 2iVD) /N. 

Using d?]) and inverting ([XTl 1X51 the time evo- 

kition can be calculated by writing \a) in terms of the 
above eigenstates and applying the time evolution oper- 
ator. Using (jAll IA3[ IA4[) again and tracing out the bath 
degrees of freedom we arrive at the reduced density ma- 
trix p(t) , which enables to evaluate the expectation value 
('^1/2 (^)) ^^'^ dynamics of the entanglement between 
the two central spins. As a measure for the entanglement 
we use the concurrence^'^ 

C(i) = max{0, VA^- VA^- VaI}, (8) 
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FIG. 5: (Color online) Entanglement dynamics for |ai) — 
and Nd = 20, resulting in pt « 6% - 30%. In the 
completely homogeneous case the amplitude is small even for 
high polarization. Generation of entanglement benefits from 
high polarization. 



FIG. 6: (Color online) Entanglement dynamics for |qi) = 
jr+) and Nd = 20, resulting in pt ~ 6% - 30%. Instead 
of an oscillating function we see discrete peaks. Variation of 
the exchange coupling has no influence because |r+) is an 
eigenstate of the central spin coupling term. 



where are the eigenvalues of the non-hermitian ma- 
trix p{t)p(t) in decreasing order. Here p{t) is given by 
{ay ® (Ty) p*{t) {(jy ® d y) , vfYieiQ p* {t) dcuotcs the com- 

plex conjugate of p{t). The coefficients c\. are of course 
products of Clebsch-Gordan coefficients, which enter the 
time evolution through the quantity 



E 



(9) 



and usually have to be calculated numerically. The main 
advantage in considering li = 1/2 is now that in this case 
a closed expression for dt can be derived^l: 



Nd\{N - Nd)KN - 2fc + 1) 
(A^ - fc + l)!fc! 



(10) 



For further details on the calculation of the time de- 
pendent reduced density matrix and the dynamical quan- 
tities derived therefrom we refer the reader to appendix 
B. 

Finally, it is a simple but remarkable difference be- 
tween our one bath system with two central spins and the 
homogeneous Gaudin model of a single central spin^^''^^, 
that even if we choose |a2) as an eigenstate and hence 
fix A; in dZl) to a single value, due to the higher number of 
eigenvalues the resulting dynamics can not be described 
by a single frequency. 



III. BASIC DYNAMICAL PROPERTIES 

We now give an overview over basic dynamical features 
of the system considered. Due to the homogeneous cou- 
plings, the dynamics of the two central spins can be read 



off from each other. Hence the following discussion of the 
dynamics will be restricted to {Sf{t)). 



A. Electron spin dynamics 

In Figs. [HO we consider the completely homogeneous 
case Jex = A/N and plot the dynamics for \a) — jfi-JJ. 
),|T+),|To) and varying polarization ph sa 2% — 30%. 
A polarization of 30% does not seem to be particularly 
high, but the behavior typical for high polarizations oc- 
curs indeed already at such a value. We omit the singlet 
case because it is an eigenstate of the system. In Fig. 
[T] the number of spins is even, whereas in Fig. [5] an 
odd number is chosen. Note that we measure the time 
t in rescaled units h/{A/2N) depending on the number 
of bath spins^. Similarly to the homogeneous Gaudin 
system-'i'"''"'"', from Figs. [TJ [5] we see that the dynamics 
for an even number of spins is periodic with a periodicity 
of TT (in rescaled time units), whereas an odd number of 
spins leads to a periodicity of 27r. This is the case for Jex 
being any integer multiple of A/N. These characteristics 
can of course be explained by analyzing the level spac- 
ings in the different situations. For example, for an even 
number of bath spins, all level spacings are even multi- 
ples of A/2N^'^; resulting in dynamics periodic with tt. 
However, if the number of spins is odd, we get even and 
odd level spacings (in units of A/2N), giving a period 
of 27T. For the given case of completely homogeneous 
couplings the dynamics can be nicely characterized: The 
number of local extrema for an even number of bath spins 
within a complete period, as well as for an odd number 
of bath spins within half a period, is in both cases given 
by — 2Ni:) + 1. This - so far empirical - rule holds for 
all initial central spin states and is illustrated in Figs. [T] 
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FIG. 7: (Color online) Spin dynamics on short time scales 
for Jex ^ 0, pb = 2/N, and |ai) = The thick solid 

lines mark the zero level (Si) = while the thick dashed 
line (lower panel) represents the threshold level (5*1) = 0.2 as 
appropriate for J ex < and small spin baths. 



andm 

Let us now investigate the spin dynamics for varying 
exchange coupling, i.e. the case Jex 7^ A/N. Note that 
for the initial central spin state \ai) = \Tq) this inhomo- 
geneity has no influence on the spin dynamics since |To) 
is an eigenstate of Si ■ S2 and 



-ffhom, Si ■ S2 



0. 



(11) 



In Fig. 13] the dynamics for \ai) — jfl-JJ-) and varying 
exchange coupling is plotted. In the upper two pan- 
els we consider the case of low polarization pt, « 10% 
for an even and an odd number of spins. The remain- 
ing two panels show the dynamics for high polarization 
Pb « 30%. In Fig. |4] the plots are ordered likewise 
for a more general linear combination of Iff-IJ.) and \To) 

, \ai) = {l/VT3) (2|M)+3|M)). 

From Figs. [31 13] we see that if the exchange coupling is 
an odd multiple of A/2N, the even-odd effect described 
above does not occur and we have periodicity of 2tt. In 
both of the aforementioned situations the time evolutions 
are symmetric with respect to the middle of the period, 
which is a consequence of the invariance of the underlying 
Hamiltonian under time reversal. For a more general 
exchange coupling, the periodicity, along with the mirror 
symmetry, of the dynamics is broken on the above time 
scales. 

Considering the case of low polarization, neither the 
dynamics of initial states with a product nor the one of 
states with an entangled central spin state dramatically 
changes if Jex is varied. However, if the polarization is 
high, the spin is oscillating with mainly one frequency 
proportional to Jex- Furthermore the amplitude of the 
oscillation is larger for the case Jex A/N than for the 
completely homogeneous case. This behaviour can be un- 



FIG. 8: (Color onhne) Position of the first zero of {Sf{t)) for 
Jex > 0, and the first intersection with the threshold level 
(^i) — 0.2 for Jex < 0, on a double logarithmic scale. We 
choose joi) = and a polarization of pt — 2/N N = 
2Nd + 2. The curves are fitted to a power law oc A''" with 
= -0.52 {Jex = (A/N)), V = -0.51 {Jex = 1.85(A/Af)), 
V = -0.53 {Jex = 0), i/ = -0.51 {Jex = -1.5(yl/iV)), v = 
-0.50 {Jex = -1.85(yl/iV)). Note that the parallel offset 
between the plots for Jex > and Jex < results from the fact 
that the intersection with the higher threshold level happens 
closer to zero. 



derstood as follows: If the polarization is high d^^-, w 1, 
whereas dfc « for A: ^ Nu. This means that calculat- 
ing the spin and entanglement dynamics, we only have 
to consider the term k = Nd. An evaluation of the co- 
effcients for the different frequencies now shows that the 
main contribution results from Etq — Es — {A/N) — Jex 
in obvious notation. Hence if the polarization is more 
and more increased, this is the only frequency left. If 
Jex = [A/N), the two associated eigenstates are degen- 
erate so that in this case the main contribution to the 
dynamics is constant. This explains why the amplitude 
of the high polarization dynamics in Figs. [31 S] is big 
compared to the one in Figs. [11 H For further details the 
reader is referred to appendix B. 



B. Entanglement dynamics 

In Figs. inilSlthe concurrence dynamics C{t) for \ai) = 
Ifl-JJ-), is plotted for the same polarizations as in Figs. 
[31 [Hand varying exchange coupling. 

It is interesting that in the second case the concur- 
rence drops to zero for certain periods of time. This is 
very similar for the case |ai) = |To) not shown above. 
As already explained concerning the spin dynamics, the 
exchange coupling Jex of course has no influence because 
|r+) is an eigenstate of Si ■ 82- 

It is an interesting fact now that for \ai) = Iff^) and 
a small polarization changing from \Jex\ > 1 to \Jex\ < 1 
increases the maximum value of the function C{t). Fur- 
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thermore we see from Fig. [5] that surprisingly the entan- 
glement is much smaller for the completely homogeneous 
case J ex = A/N than for J^x 7^ A/N even for low polar- 
ization. 



IV. DECOHERENCE AND ITS 
QUANTIFICATION 

Depending on the choice of the exchange coupling, the 
dynamics of the one bath model can either be symmetric 
and periodic or without any regularities. It is now not 
entirely obvious to determine in how far these dynamics 
constitute a process of decoherence. Considering for ex- 
ample the spin dynamics for an integer Jex and an even 
number of bath spins shown in Fig. [1] one can either 
regard the decay of the spin as decoherence or, especially 
due to the symmetry of the function, as part of a simple 
periodic motion. In Ref.^^ the first zero of {Sf{t)) has 
been considered as a measure for the decoherence time. 
In Fig. [7] we illustrate examples of the spin dynamics 
on short time scales for Jex > 0, Jex < and a vary- 
ing number of bath spins. For J^x > this procedure 
is straightforward meaning that {Sf{t)) crosses the hor- 
izontal line (Sf) = before reaching its first minimum 
with {Sf{t)) < 0. However, for Jex < and a sufficiently 
small number of bath spins, as seen from the lower panel 
of Fig. [71 such a first minimum is attained before the first 
actual zero (5f (i)) — 0. This first zero occurs indeed at 
much large times t whose scaling behavior as a function 
of system size N is clearly different from the zero posi- 
tions found for J^x > 0, as we have checked in a detailed 
analysis. Thus, our evaluation scheme needs to be mod- 
ified for Jex < 0. An obvious way out of this problem 
is to either consider large enough spin baths where such 
an effect does not occur, or to evaluate the intersection 
with alternative "threshold level" (5^ ) > 0. In Fig. [7]we 
have chosen (Sf) = 0.2, which will be the basis of our fol- 
lowing investigation. As a further alternative, one could 
also consider the position of the first minimum of {Sf{t)). 
Hence, strictly speaking, it is not per se the first zero of 
{Sf{t)) < which is a measure for the decoherence time, 
but the scaling behavior of the dynamics on short time 
scales. Following the route described above, in Fig. [5] we 
plot the positions (measured in units of h/ {A/2N)) of the 
first zeroes of {Sf{t)) for Jex > 0, and of the first intersec- 
tions with the threshold level shown in Fig. [7]for Jex < 0, 
on a double logarithmic scale. We choose a weakly po- 
larized bath N — 2No + 2 ^ pi, = 2/N, approaching 
the completely unpolarized case for N ^ 00. The ab- 
solute values of the positions for Jex > and Jex < 
differ slightly from each other, which results from the 
fact that the intersection with the threshold level at 0.2 
happens closer to zero than with the usual threshold level 
(Si) = 0. Nevertheless, the scaling behavior is very sim- 
ilar in all cases, and each curve can nicely be fitted by a 
power law oc {N + 2Y with v k, —0.5, a result similar to 
the one found for the homogeneous Gaudin system with 



only one central spin--. 

In a GaAs quantum dot the electron spins usually in- 
teract with approximately = 10^ nuclei. Assuming 
the hyperfine coupling strength to be of the order of 
A = 10~^eV, as realistic for GaAs quantum dots^^, this 
results in a time scale of Nh/{nA) ~ 1.31 • lO^'^s. If we 
now use the above scaling behaviour + 2, we get 

a decoherence time of 131ns, which fits quite well with 
the experimental datap^iii^i^. This is an interesting re- 
sult not only with respect to the validity of our model: 
As explained following equation ©, generally decoher- 
ence results "directly" from the electron-nuclear flip-fiop 
terms and through the superposition of product states 
from the z terms. Above we calculate the decoherence 
time for \ai) = Iff-IJ.), where the infiuence of the z terms 
is eliminated. The fact that we are able to reproduce 
the decoherence times suggests that the decoherence time 
caused by the fiip-fiop terms is equal or smaller than the 
one resulting from the z parts of the Hamiltonian. It 
should be stressed that we calculate the decoherence time 
of an individual electron T2 here. In Ref<^ the decoher- 
ence time of an ensemble of dots Tj* been calculated 
yielding Ins for a GaAs quantum dot with 10^ nuclear 
spins. 

It is now a well-known fact for the Gaudin system that 
the decaying part of the dynamics decreases with increas- 
ing polarization '^l. A numerical evaluation shows that 
this is also the case for two central spins. As explained 
in the context of Figs. [TJ [21 O S] the oscillations of our 
one bath model become more and more coherent with 
increasing polarization. Together with the above results 
for the decoherence this means that, although the ho- 
mogeneous couplings are a strong simplification of the 
physical reality, our homogeneous coupling model shows 
rather realistic dynamical characteristics on the relevent 
time scales. This is plausible because artifacts of the ho- 
mogeneous couplings, like the periodic revivals, set in on 
longer time scales. 



V. CONCLUSION 

In conclusion we have studied in detail the hyperfine in- 
duced spin and entanglement dynamics of a model with 
homogeneous hyperfine coupling constants and varying 
exchange coupling, based on an exact analytical calcula- 
tion. 

We found the dynamics to be periodic and symmetric 
for Jex being an integer multiple of A/N or an odd multi- 
ple of A/2N , where the period depents on the number of 
bath spins. We explained this periodicity by analyzing 
the level spectrum. For Jex = A/N we found an em- 
pirical rule which charaterizes the dynamics for varying 
polarization. We have seen that for low polarizations the 
exchange coupling has no significant influence, whereas 
in the high polarization case the dynamics mainly con- 
sists of one single frequency proportional to Jex ■ It is not 
possible to entangle the central spins completely in the 
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setup considered in this article. 

Following Refj^ we extracted the decoherence time 
by analyzing the scaling behaviour of the first zero. In 
the case of negative exchange coupling the dynamics 
strongly changes on short time scales and instead of 
the first zero we considered the intersection of the 
dynamics with another threshold level parallel to the 
time axis. Both cases yield the same result which is 
in good agreement with experimental data. Hence 
the scaling behaviour of the short time dynamics can 
be regarded as a good indicator for the decoherence time. 



As already mentioned in the text, the states are labelled 
by the quantum numbers J, m, S corresponding to the 
operators , J^^ , . The rest of the quantum numbers 
due to a certain Clebsch-Gordan decomposition of the 
bath is omitted. For the eigenstates of the central spin 
term Si ■ S2 we used the standard notation: 
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Appendix A: Diagonalization of the homogeneous 
coupling model 

The eigenstates of -ffhom can be found directly by iter- 
ating the well known expressions'^^ for coupling an arbi- 
trary spin to a spin S = 1/2. Two of these states lie in 
the triplet sector: 



\T+) 
\To) 
\T-) 
\S) 



m) 

^(IM)-IM)) 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



I + m + 1 I + m , 
|/+l,m,l) = J———.-——\I,rn-l)\T+) 



21 + 2 21 + 1' 



I + m+l I -m+l 



I + l 



21+1 



I,m)\To) 



7 — 771+1 I — m , 



21 + 2 21+1 



\Lm + l)\T) 



(Ala) 



\I-^,n.^) = ^r-^■'-^:;^\I,rn-lm 



21 



21+1 



I — m I + m 



21 +V 



Lm)\T, 







ll + m 1 + 777-1-1 



21 



21+1 



7, 777 + i)|r_ 



(Alb) 



The remaining two eigenstates are superpositions of sin- 
glet and triplet states. As the expressions are rather 
cumbersome, it is convenient to introduce the following 
notation in order to abbreviate the Clebsch-Gordan co- 
efficients: 



r 2 2 2 2^ 



I + m 7 — 777 + 1 I + m I + m I — m I — m I — m 7 + 777 + 1 



27 



27 + 1 



27 27 + 1 ' V 27 27 + 1 



27 



27+1 



7—777+1 I + m 

21 + 2 ' 27+1' 



1 7 



27 + 2 



27+1 



1 7 + 777 + 1 /7 + 777 + 1 7 - 777 



27 + 2 



27+1 



27 + 2 27+1 



With this definitions the superposition states can be 
written as: 
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|1) = ^J|/,™-l)|T+) + ^^i-^|/,m)|ro> 
- ^i|/,™+l)|T_) + ^^i^|/,™)|5) 



|2) = /if|/,™-l)|r+) + 



72 



Aii|/,m+l)|T_)- 



These states are degenerate with respect to 77hom, hence 
we are left with the simple task to find a superposition 
of |1) and |2), which eliminates \I,m)\S). Obviously this 
is given by 

where Nt = a/— (/ + 1)^^ + + 4 is the normaliza- 
tion constant. Inserting |1) and |2) this reads: 



|/,m,l) 



-1)|T+ 



/,m)|rr 



0; 



Together with the singlet state 

|/,m,0) = |/,m)|5) 



(A3) 



(A4) 



this solves our problem of diagonalizing the one bath ho- 
mogeneous coupling Hamiltonian. Furthermore (jAip and 
(jA3p give a solution to the very general problem of cou- 
pling an arbitrary spin to a spin S ^ 1. 



The time evolution of the initial state results from the 
application of the time evolution operator U = e~'K^*. 
It follows: 

\a{t)){a{t)\ = \Ua){Ua\ 

= y^Q!,o;*|[/-0,)([/-0j-| 



\^,){^,\ (B2) 



As Oi acts on Hi, the other degrees of freedom have to 
be traced out 

p,(i)-Tr„\„,(|a(t))(a(i)|), 
finally giving the time evolution of the operator: 



(0,(i)) =Tr«, (A(t)O,) 



(B3) 



Usually such calculations are done numerically, but for 
our homogeneous coupling model it is possible to derive 
exact analytical expressions for the dynamics of the two 
central spins. 

Following the general scheme, we have to write the 
initial state in terms of energy eigenstates first. As ex- 
plained in the text, we consider |a) = |ai)|Q;2), where 
\ai) is an arbitrary central spin state and \a2) is a prod- 
uct state in the bath Hilbertspace Hn- Using ([7]) it fol- 
lows: 

\m)\a,} = 2 ^ 'Vi>iy - f - Nn,{S^}) (B4) 

fc=0 {Si} 

The eigenstates (|A11 IA31 IA4p are given in terms of prod- 
uct states between a basis element from (jA2l) and an 
eigenstate. Hence we can find the coefficients of (IBip by 
solving (jAl[ IA31 IA4p for these states and inserting them 
into (|B4I) . If we arrange the coefiicients from (jAll IA31 
4| into a 4 X 4 matrix V according to 



Appendix B: Calculation of the time-dependent 
reduced density matrix 

Let be a time-independent Hamiltonian acting on 
a product Hilbert space V. — ^f^i'Hi- We denote its 
eigenvectors by ji/'i) and the corresponding eigenvalues 
by Ei. In the following we calculate the time-dependent 
reduced density matrix for an initial state which is a pure 
state and derive the time evolution {Oi{t)) associated 
with an operator Oi acting on T-Li. Then we consider 
the Hamiltonian ([3]) and give some more details on the 
corresponding calculations for our model. 

As the eigenstates of H span the whole Hilbertspace 
the initial state |a) of the system described by H can 
be written as 

|a)=^a,|V^). (Bl) 

i 





in) in \T'-) 


\S')\ 


|/ + l,m,l) 






|/,TO,1) 






|/-l,m,l) 






|/,m,0) 




■■• / 



this is simply done by transposing V. Here \T'j^) — 
|/ — 1, to) |r+) and analogously for the other states. In or- 
der to abbreviate the following expressions we denote the 
energy eigenstates by \i}}i) as in the general considerations 
above and number with respect to (IB5|) . Analogously we 
introduce the shorthand notation \i) for the basis states 

(EH. 

In order to avoid further coefficients we choose \ai) 
to be the j-th element of (IA2[) and find the following 
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expression for the decomposition of the initial state into 
energy eigenstates 

4 Nd 

b)i"2) = EEE4'''^^TiV'z), (B6) 

1=1 k=0 {SJ 

where it is has to be noted that the elements and 
the eigenstates \ipi) depent on the quantum numbers the 
sums run over. Hence in our case the coefficients ai and 
the eigenstates in fact have more than one index. 
Inserting ((B6)) and (|XT1 \Mi EH) in ([B2|) and tracing 
out the bath degrees of freedom, we finally arrive at the 
reduced density matrix of the two central spins 

p(0:=Tr«„(|a(i))(a(t)|)= (B7) 

fe=0{Si} l,m,n,o=l 
dk 



If we now choose Oi = Sf, we have to trace out the 
second central spin. Inserting the result into (|B3|) then 
gives rise to the time evolution {Sf{t)). This is given 
by (|B7|) with n — o, multiplied by coefficients resulting 
from the eigenvalues of 5f . As mentioned in the text, 
for high polarizations « if fc 7^ Njj. Fixing l,m we 
can calculate the contribution of the respective frequency 
by evaluating the remaining sum over n. If the polariza- 
tion is strongly increased, all frequencies are suppressed 
except for E2 — E4 = Etq — Es ■ 
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